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Several technical aspects concerning the evaluation of the overlap between two mean field wave 
functions of the Hartree Fock Bogoliubov type, are discussed. The limit when several orbitals 
become fully occupied is derived as well as the formula to reduce the dimensionality of the problem 
when exactly empty orbitals are present. The formalism is also extended to deal with the case 
where the bases of each of the wave functions are different. Several practical results concerning the 
evaluation of pfaffians as well as the canonical decomposition of norm overlaps are also discussed in 
the appendices. 
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I. INTRODUCTION 

In a recent publication [1] a new formulation, based on 
the pfafHan of a skew-symmetric matrix, has been pro- 
posed to compute the overlap between mean field wave 
functions of the Hartree Fock Bogoliubov (HFB) type, in- 
cluding its phase (or sign, for real overlaps). The result 
for the overlap is obtained by recursing to the powerful 
concept of fermion coherent states 0-Q and it involves 
the evaluation of a quantity called the pfaffian of a skew- 
symmetric matrix -see. for instance [7] for a definition of 
the pfaffian in a physical context- that is similar in spirit 
(linear combinations of products of matrix elements) to 
the determinant of a general matrix. In the derivation of 
the formula for the overlap (H it is assumed that the two 
mean field wave functions (of the HFB type -see ^ for 
definition and properties-) can be related to a common 
reference one (usually chosen as the true particle vacuum) 
by means of their Thouless parametrization. It is very 
likely to find cases where the Thouless parametrization 
is ill defined because, for that particular case, it involves 
the inverse of a (near) singular matrix. This situation 
corresponds to the presence of particles (or quasiparti- 
cles in the general case) that have an occupancy of one, 
rendering the wave function orthogonal to the common 
(or reference) wave function. To handle those singular 
cases it was suggested in [l| to just change the common 
reference mean field wave function in order to modify 
the occupancies with the hope that none of them will 
be close to one. However, it is desirable to have an al- 
ternative for those cases where the change of reference 
wave function is either not possible or too cumbersome 
to carry out. Therefore, I have considered the formal 
limit of occupancies going to one and I have obtained 
a formula which is well defined in that limit and pro- 
vides a sound answer to such singular limit. By using 
the same kind of ideas I have also handled explicitly the 
case where some particles (or quasiparticles) have zero 
occupancy and therefore they do not contribute to the 
overlap. The formula obtained in this case involves ma- 
trices of smaller dimensions than the original ones and 



' lluis.robledo@uam.esl 



therefore should be regarded as a thrifty alternative to 
the original formula for situations where the wave func- 
tions are expanded in huge bases. Similar manipulations 
to reduce the size of matrices and/or deal with full occu- 
pancies, have been considered in Refs 0-[l2|. I also ad- 
dress the case where each of the mean field wave functions 
are expressed in different single particle basis related to 
each other by a general (not necessarily unitary) trans- 
formation. The result obtained is useful, for instance, to 
compute the overlap of the operators for spatial transfor- 
mations (as translations or rotations) between arbitrary 
mean field HFB wave functions. The result is general 
enough as to allow for transformations that do not map 
the single particle basis into itself (non complete basis 
under the transformation). This was already considered 
in [9] in a general framework but not considering the 
present formulation including the phase of the overlap. 
In [12] the implications of considering two different bases 
are also address, but there it is implicitly assumed that 
both bases share the same block structure defining the 
conjugate states. This is a limiting assumption that do 
not hold in general for time reversal violating (crank- 
ing, for instance) wave functions. Finally, some useful 
results concerning the formal evaluation of pfaffians are 
discussed in the appendices. These results can be of inter- 
est in other branches of physics where the use of pfaffians 
is becoming increasingly popular (Tsl - fTsi . To cover also 
more practical aspects, the reader is referred to [3 for a 
thorough description of useful algorithms to compute nu- 
merically and symbolically the pfaffian of arbitrary skew 
symmetric matrices. 

The relevance of the results present here is a direct 
consequence of the increasing popularity of the so called 
"beyond mean field methods" in nuclear physics (Tol . [itI - 
[iH that demand the evaluation of both the modulus and 
phase of the overlaps between arbitrary HFB wave func- 
tions. A reliable determination of the sign of the norm 
can also be useful in to order to pin down the location 
of the zeros of the HFB overlaps [l^l- This determina- 
tion would eventually be useful to get rid of the so called 
"pole problem" that plagues present beyond mean field 
calculations [23 - [25| . 

In section |TT] the formulas pertaining the two limits 
considered are derived and their implications discussed. 
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We also explicitly show how to implement the change 
to a different common reference HFB state that could 
be an easy alternative in some cases. In section llllt he 
case where the two HFB wave functions are referred to 
different bases is discussed. Finally, in appendices IIIIl HIl 
ancJl]some relevant results required in the derivations are 
discussed. 



II. EVALUATION OF THE OVERLAP IN 
VARIOUS LIMITS 

To begin with, it is convenient to recall the nota- 
tion used in [l| as well as the formulas obtained there. 
The goal is to evaluate both the modulus and phase of 
the overlap {(j)o\4>i) between two HFB wave functions 
\(f)()) and \(f)i) which are given in terms of the Thouless 
parametrization of a general HFB wave function 
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The diagonal matrix M^^^has matrix elements 
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The extreme values of the ratios Vj /u^ are infinity for 
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fully occupied levels (w^'^ = 1) or zero for empty levels 
(ii(') = 0). For further developments it is convenient to 
single out those values and write 



Ni^ 

o(*: 



(7) 



where the diagonal matrix O*^*' contains the if(') diagonal 
elements belonging to the extreme values, infinity or zero, 
mentioned above. The dimension of this matrix is K^^^ x 
if*^'). Each of the two limiting cases require different 
considerations and hence we will from now on considered 
them separately. 



A. Limit of fully occupied levels 



The skew-symmetric matrices M^*) = of 
dimension N x N (N is assumed to be an even number 
N — 2q a.s required for fermions), are related to the co- 
efficients U^'^^ and of the Bogoliubov transformation 
defining the quasiparticle annihilation operators 



associated to \(j>i). The above wave functions are not 
normalized to one, but as to have {0\(j)i) = 1 instead. As 
shown in [l| the overlap can be written as 



with 



/,o|0i) = (-l)^(^+i)/2pfM 
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is a 2N X 2N matrix. To obtain the results of the present 
section the Bloch-Messiah decomposition of the Bogoli- 
ubov amplitudes [g, namely = C^*) and 
= £)(*)*\/(0(7(«)^ ig used. In the previous expres- 
sions, and C(') are given unitary matrices and C/'^*-' 
and V^'^ are real matrices with special diagonal forms. 
By using this decomposition we can write 



(4) 
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where the skew-symmetric matrix is in "skew- 

symmetric diagonal" (or canonical) form 
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In this case, there are K^^^ fully occupied levels in each 
of the HFB wave functions and the corresponding 

diagonal elements of the matrices M (the ones cor- 
responding to O'*-' in Eq. ([7])) tend to infinity. This 
is a serious challenge, as the overlap Eq ([2]) as well as 
the norm of the diverge. The divergence has to 
be regularized and singled out of the overlap in order 
to cancel it out with the diverging factors coming from 
the norms of the HFB wave functions. To this end we 
write M, 

the "canonical regularized" (CR) matrix 
/ 



= R<-^^ M^Ir<-^^^ where we have introduced 
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as well as 
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In all the cases I^(i) represents the unit matrix of dimen- 
sion K'^^'^ . With the above definitions the matrix M of Eq 
([3]), which enters the expression of Eq © for the overlap, 
is factorized as 
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and 
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Using now the property pi{B^AB) — dct(i?)pf(A) we 
obtain 

pf(M) = det(i:'(i))det(i:'(")*)det(i?(i))det(i?(")*)pf(M) 
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In the case of fully occupied levels, the diagonal 
matricesO'-*-' introduced in Eq. ([5]) tend to infinity. As 
a consequence, the determinants det(i?'^*^) diverge as do 
some of the matrix elements of w'^\ The later is not a 
problem as just the inverse of this matrix is required in 
Eq. ® 

and -1 is a well defined quantity 



As a consequence of this structure, the matrix S of Eq 
^ is, in this limit, a matrix where the last K^^^ rows 
and last isT^"^ columns are set strictly to zero. We can 
use this property together with the special structure of 

(i) 

the matrices M^]^ to reduce the size of the matrices to 
be considered in the evaluation of the pfafhan. This will 
be addressed in the next subsections in a slightly differ- 
ent context. The only truly diverging quantities, namely 
the determinants det(i?'*-'), cancel out when we compute 

the normahzed overlap (</'oki)/\/|(<<5o|<Po)PI(vi|'/'i)P 
as, from the previous formulas, is proportional 

to |det(i?W)|2. 



occupancy. This is an advantage as the number of non- 
zero occupancy levels is usually much smaller than the 
total dimensionality of the basis used. I will use the same 
notation as before, but in this case, the O(^) are diago- 
nal matrices whose diagonal matrix elements are made 
to tend to zero. The number of diagonal elements, K'-'^\ 
represent the number of empty levels in each HFB wave 
function \4>i). In this case it is convenient to reorder the 
single particle basis by using unitary similarity transfor- 
mations P23 (see appendix H)) permuting blocks 2 and 3, 
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where the matrix N^^^ has dimension 2{N — K^^^)x 2{N~ 
K^'^) and O^'^ has dimension 2K^'^ x 2K^'\ The matrix 

(i) 

carrying out the permutation depends on the basis 
as the size of each block TV*^*^ is not necessarily the same. 
The impact of this dependence will show up in the eval- 

(i) 

nation of the pfaffian where the determinant of is 



required. Now, the new unitary matrix D 



'CR^R 
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Now a technical detail related to the way how to handle 
the situation where the bipartite structure of M^*'^ is dif- 
ferent from the one of M^]^ has to be considered. The dif- 
ferent bipartite structure happens when the dimensions 
and if^i) differ, or in other words, that the number 
of empty levels in each mean field wave function is differ- 
ent. The optimal situation in terms of simplicity is that 
when both JsT^^^ are the same. This is the strategy we 
will use in the following by assuming the same bipartite 
structure in both M^^^ and D^^ with a common dimen- 
sion Ks which is the smaller of ivT^") and isT'^^). With 
this in mind, we endow the matrices M^']j, with the new 
bipartite structure 



B. Limit of fully empty levels 

Another situation often encounter in numerical ap- 
plications is when many orbitals have zero occupancies 

— and therefore their contribution to the overlap is 
zero. To disentangle those contributions and reduce, in 
this way, the computational cost of the evaluation of the 
norm it is convenient to consider the limit of fully empty 
levels. This limit has been considered by other authors 
[^-[1^ in the past but in the more traditional formula- 
tion of the overlap not including the phase. In the limit 
of fully empty levels I will show that the evaluation of 
the overlap can be reduced to considering matrices with 
dimension equal to the number of orbitals with non-zero 




(11) 



where n'^'^ me 2{N - Ks) x 2{N - Ks) matrices and the 
0(*) have the same dimension 2Ks x 2Ks- In the same 
way, we will use below the same bipartite structure for 



the matrix D 



R 



^12 



(12) 



22 



where D^^ is a square matrix of dimension 2 (TV — Ks) x 
2{N-Ks), D^^ is of dimension 2{N-Ks) x 2Ks , d!^^ 
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is of dimension 2Ks x 2(A^ — Ks) and finally is of 

dimension 2Ks x 2Ks- If we now use Eqs. pT|) and 

to reconstruct the matrix M*^'^ of Eq. (|10p we realize 

that in the O^'^ — >■ limit the sub-matrices ^^i^ -^22* 
do not enter the final expression of M^'^\ We can use 
this freedom to chose those "arbitrary" matrices in such 
a way as to simplify some of the results to be obtained 
below. A possible, and convenient, choice is D^l = and 



D 



12^3 • The matrix obtained with this choice will 



(0 



be denoted by Dj^ with 



D 



D 







11 

D21 hKs 



(13) 



and has the nice property of having a simple inverse 



D 



D^^}-' 



involving only the inverse of the matrices D^^ which have 
a moderate dimensionality. This is the property guiding 
the choice made, as we are implicitly assuming that N -— 
Ks ^ Ks and it has to be kept in mind that the cost of 
most of the matrix operations grow as the cubic power 
of their dimension. 

With all these definitions the matrix M which enters 
the formula for the overlap is written as. 
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where 
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^(l)-l /^(O) 
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matri- 



(please note that the unitary character of the D 

ces is lost with the introduction of the Dj^ ones). Using 
now the properties of the pfafhan we obtain 



pf(M) = det(L)^^Vet(-C)2^*)pf(M) 
= det(Z?ilVet(i^ii^*)pf(M) 



with 
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Let us now analyze the structure of the block matrix U 
entering the definition of M. Using Eq p4| together with 
Eq (fT3| allows to write 



U 



Un U12 
U21 U22 



with 



t/n=(<Mli)" 

t/l2--C/llM°) + 



U21 
U22 



'R 21 
-^i? 21^11 
+ ^R21^^l^R21 



(18) 

(19) 
(20) 
(21) 



With this definition, the matrix M of Eq ([28| acquires, 
in the limit where O''' — 0, the block structure 








Un 


U12 








U21 


U22 


-c/fi 


-ul 







-uL 


-ul 









This block structure is still not beneficial for the simpli- 
fication of the corresponding pfaffian and we need to use 
the exchange matrices defined in appendix U to exchange 
blocks 2 and 3 



= P23MPI^ = 



Using the formulas of appendix |TT] for the pfaffian of a 
bipartite matrix we obtain 



-f/fl 


Un 




~ul 


U12 






.-UL 


C/21 





~uL 


c/22 





Pf(Mfl.) = Pf 



X Pf 




where 



^12=C/ll+C/l2C/2;lC/21 



(22) 



The first pfaffian in the right hand side of the above ex- 
pression is simply given by (—1)^^ det C/22 whereas the 
second pfaffian can be computed using again the expres- 
sion for the pfaffian of a bipartite matrix. Collecting all 
the terms together we finally obtain 



Pf(MK) - (-1) 



Ks 



X Pf -N. 



det(C/22)Pf 

(0) 



"c 



r(l)-l 



X^N^ ' X12 



Taking into account that det(P23) = 1 we finally obtain 
pf(M) = (~1)^- det{D[\^)det{D^°^^ *) det(C/22)Pf [n^c^) 

'X12) (23) 



X Pf (-iV^°)' 



Xr2A^«- 



The advantage of this ugly result over the general ex- 

^11 



pression is that the dimensionality of the NX. , iV, 



(1) T^m 



and Un matrices is 2(A^ — Ks) which is much smaller 
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than the one of the original problem (27V). The only big 
matrix in Eq. (j22p is the inverse of U22 with dimension 
2Ks, however, this does not pose a challenge as its spe- 
cial structure Eq. (PT|) is very well adapted to the use of 
the Sherman-Morrison formulas for the determinant and 
inverse of this kind of special matrices [l^ . 



C. Using a different reference vacuum 

From the above discussion, it is clear that the structure 
and properties of the M^'''> matrices is intimately related 
to the reference vacuum used to express the HFB wave 
functions This suggests to use another reference vac- 
uum, instead of the true vacuum implicitly assumed in 
the previous discussions, with the hope that the new ma- 

trices M will acquire a structure where there would be 
no fully occupied quasiparticles and the number of empty 
ones will be very large and comparable to the size of the 
basis. Matrices with that properties will not require the 
use of the "fully occupied limit" formulas and will benefit 
from the reduction in computational burden of the "fully 
empty limit" results. To be more specific, let us consider 
the example where the \<j)i) correspond to two HFB wave 
functions with the same average of the number of par- 
ticles and different quadrupole deformation parameters 

(i) 

^2 • It looks rather intuitive that a more appropriate ref- 
erence vacuum than the true vacuum could be the HFB 
wave function |0) with the same average of the number 
of particles and a deformation parameter 52 close to both 
(for instance the mean value ^{q^'^ +92^''))- It is to be 
expected that, with respect to this reference vacuum 

the new U amplitudes will be very close to the identity 

(i) 

matrix whereas the new V amplitudes will be small. 
In other words, the HFB wave functions will be rep- 
resented by a linear combination of quasiparticle exci- 
tations of the reference state \(f)) with small amplitudes 
that quickly decrease with the number of quasiparticle 
excitations (i.e. the amplitudes of the four quasiparticle 
excitations much smaller that the amplitudes of the two 
quasiparticle ones). The expected properties of the new 

U and V amplitudes imply that most of the quasi- 
particle excitations will correspond to the "fully empty" 
limit discussed in the previous subsection (with the as- 
sociated advantages of having to deal with matrices of 
small size) and far from the problematic fully occupied 
limit of subsection III Al 

Let us consider the new reference vacuum |0) with the 
associated creation Q;~'"and annihilation a quasiparticle 
operators which are defined in terms of a single particle 
basis of creation and annihilation operators by means of 
linear combinations involving the U and V amplitudes 



a 



a 
+ 




ators a^'^ and a^*^"*" with amplitudes W^^\ Using the 
unitarity of the matrices T4^(*) and W we can express the 
set of quasiparticle operators a^*^ and aW+ in terms of 
the d and d'^ ones as 



(0 



a 



a 




Using Thouless theorem we can also express the wave 
functions ( satisfying (^|0('^) = 1 and therefore different 
from the previous \(j)i) by a normalization factor) in terms 
of the 1 6) reference vacuum 




(24) 



with Af(*) = (y(0[/(i)-i)* (not to be confused with the 
diagonal matrix of Eq (P5)l ). These two matrices can be 
easily computed once the W coefficients have been given 
and it is even possible to give an analytical expression 
[27! in terms of M and M^*) 



M(*) = Q(M^ -hMW)(I + M+MW)-l(Q+)- 



(25) 



with Q={1 + MM+)-i/2. AH the formulas given above 
(and below) are equally valid for the wave functions given 
in the form of Eq p4)) with the amplitudes of Eq. ([25)) . 



The same relation holds true for the quasiparticle oper- 



III. DIFFERENT SINGLE PARTICLE BASES 

It is very common that the HFB wave functions 
are defined in terms of different single particle basis, with 
creation and annihilation operators a^(«) and ak{i) that 
will carry indexes [i) to indicate the HFB wave function 
they belong to. Those bases are usually not complete and 
therefore they span different subspaces of the full Hilbert 
space. As a consequence, the formulas obtained above 
cannot be used because they implicitly rely on a common 
basis for the two HFB wave functions [9]. The strategy 
to overcome this problem is to find a bigger subspace en- 
compassing both subspace and use an orthogonal basis 
defined in that bigger subspace for the two HFB wave 
functions In an early consideration of this prob- 

lem [9] I used the whole Hilbert space as the common 
subspace. Another possibility, explored in this paper, is 
to consider the subspace generated by the union of the 
two subspaces. In this case, special care has to be taken 
with the resulting basis, union of the of the two original 
bases, as it can be redundant (i.e. it can contain linearly 
dependent vectors). 

Let me start by considering the two basis {a^{0), k ~ 
0,...,7V(0)} and {o+(l), k = 0, ...,7V(i)} which are de- 
fined in terms of single particle creation operators (typi- 
cally those of the harmonic oscillator basis) and with di- 
mensions iV(o) and ^^(1), respectively, which do not neces- 
sarily coincide. It is worth introducing the set of creation 
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operators 



At 



a+(0) /I = fc A: = l,...,7V(o) 



embracing the two sets of creation operators of the bases. 
They satisfy the commutation relations {A^, A^} — M^u 
and {Ap,, A,j} = A^} = where the overlap matrix 

A/" (dimension (A^(o) + ^(i)) x (^(o) + ^(i)) is given in 
terms of the rectangular matrix Tk y — {afc(0)5 o-^i (1)} — 
(— |afc(0)a^, (1)|— ) by the expression 



(26) 



The overlap matrix is hermitian, semi-positive definite 
and therefore can be diagonalized by a unitary transfor- 
mation £), i.e. M = -DnD+where the diagonal matrix n 
of the eigenvalues is of dimension (iV^o) -|- N^i-^) x {N^q) + 
iV(i)). In order to deal with the zero (or smaller than a 
given threshold) eigenvalues case (which signals the ap- 
pearance of linearly dependent basis states) it is conve- 
nient to introduce the notation 





where n is a diagonal matrix with the eigenvalues dif- 
ferent from zero and e is the diagonal matrix of dimen- 
sion containing those eigenvalues with value zero (or 
smaller than a numerical threshold). It is convenient in 
the ensuing developments to consider that the matrix e 
is different from zero and therefore can be inverted. At 
the end of the calculations e will be made to tend to 
zero to obtain the final result. Taking this regulariza- 
tion scheme into account, we can define the square root 
of the overlap matrix A/"^/^ = Dn^/^D+ and its inverse 
M~^l'^ = DnT^/'^D^ that are required to define the op- 
erators 

B, = Y.N;:y^A, 

V 

as well as the inverse relation = Yliv^tJ-f ^t^- 

The creation and annihilation operators i?^ and B^, 
satisfy canonical commutation relations {_B^,_B+} — 



-1/21 



They are introduced to ex- 



press the HFB wave functions of Eq. (|24)) in the standard 
way as 

|0,)=expJ^iiV«B+i?+ i|0) 



with the matrices of dimension (A^(o) + -^(i)) x (^(o) 
^(1)) 



2 * 



(27) 



given in terms of the extended matrices (also of dimen- 
sion (Af(o) + iV(i)) X (iV(o) + iV(i))) 






Af(i) 



(28) 



As the operators i?+ and satisfy canonical commu- 
tation relations, it is now possible to apply the standard 
formalism already developed in Ref [l| to write 

= ^W(o)+W(i,pf(M) 

where the matrix M entering the argument of the pfaffian 
is given in terms of the matrices defined in Eq ((28)) as 




Using the results of aDpendix lIIII concerning the eigenval- 
ues and eigenvectors of the norm matrix J\f (to simplify 
the notation in the following we consider the size of the 
two basis to be equal iV(o) — iV(i) = N) and using the 
definition of Eq (P7)) we have 



n 



1/2 



n 





1/2 



which defines the matrices N^j^^ as 



N 



D 



4' 












1/2 



Please note that these matrices are well defined when 
some of the eigenvalues of the norm overlap matrix M go 
to zero (i.e. some of the elements of the diagonal matrix 
n. are zero). Using the explicit form of the matrix D 
given in appendix IIIII in terms of the matrices E and F 
entering the Singular Value Decomposition (SVD) of the 



matrix T and defining the auxiliary matrices E± 

1/2 

and F± = F we get 



1/2 



E 



(0) 



N 



and 




M(o) _M(o) \ / El 







N 



(1) 




E^ 




Finally, using known properties of the pfaffian we can 
express the overlap in terms of the iV]^ as 



(1) 



-N 



(0). 
D 



(29) 



7 



which is the final expression for the overlap. 

To finish this section it is worth considering the limit 
where the two bases are connected by means of an unitary 
transformation (see appendix IIIip . In this case, the SVD 

1/2 

of T is trivial and we have 
and F — T. Using these values we have 



V2I, n^J^ ^0, E = 



paid to the problem of redundancy of the enlarged sub- 
space. The formulas given in this paper are a practical 
complement of the general one given in [l| and should be 
useful for a practical implementation of the calculation 
of the overlaps of two different HFB wave functions in 
the most general case. 



^(1) ^ / rAf(i)T^ 
^ I 
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As a consequence of the zero eigenvalues of the overlap 
matrix, the matrices acquire a bipartite structure 
where only the upper diagonal block is different from 
zero. In this case, the matrix in the pfaffian in Eq. 
becomes a block matrix 



/TM(i)r^ 









V 











\ 

-I 








/ 



The pfaffian of this matrix can be simplified by using 
the results of appendix U and, after exchanging blocks 2 
and 3, bring the matrix to block diagonal form. Once in 
block diagonal form, the pfaffian can be reduced to the 
product of the pfaffian of each of the diagonal blocks by 
using a simplified version of the results of appendix HIl for 
the pfaffian of a bipartite matrix. The final results is 



_l)A^(W+l)/2pf j 

as expected (see Ref 
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I. REORDERING OF MATRICES AND ITS 
IMPACT IN THE PFAFFIAN 

In some situations we will have to make use of eventual 
peculiarities of the block structure of the skew-symmetric 
matrix in order to simplify the final expression of the 
pfaffian. To do so, it is convenient to know how to reorder 
rows and columns of a matrix as well as the impact of 
such reordering in the pfaffian. A useful set of matrices 
is the one of the matrices E{i,i) that exchanges columns 
i and j of a matrix when multiplied to the right hand 
side of that matrix, i.e. 



-1 o,j a^+i 



The matrices£'(ij) are characterized by the matrix ele- 
ments 



IV. CONCLUSIONS 

We have analyzed the formula for the evaluation of the 
norm overlap of two different HFB wave functions in the 
situation where some of the occupancies of the quasipar- 
ticle levels are one and the standard approach leads to 
indeterminacies what have to be singled out in order to 
obtain a meaningful answer. In the case where there are 
fully empty single particle levels the overlap formula is 
well behaved but this limit is also addressed as a (sig- 
nificant) reduction in the computational burden can be 
obtained if the number of empty levels is large enough. 
Finally, the case where each of the two HFB wave func- 
tions are expressed in different single particle basis is ad- 
dressed and the formalism to compute the overlap in this 
situation is developed. A common basis, union of the 
other two is considered, and special attention has to be 



E{i,j)ki = Ski - SkiSii - SkjSij + SkjSii + SkiSij 

and are unit matrices where the elements i and j of the 
diagonal are set to zero and the elements i, j and j, i are 
set to one. 

Another useful set of matrices is the one of the S{i,j) 
matrices such that, when multiplied to the right hand 
side of a matrix, moves the column j of the matrix to the 
position of column i {i < j) and then shifts column i to 
position i + 1, column i -|- 1 to position i -|- 2, and so on, 
up to column j — 1 that is shifted to column j', i.e. 
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In terms of matrix elements they are given by 

s=i s=i 

These matrices are unit matrices where the 1 in position 
i,i is shifted to position , the 1 in position 

is shifted to position i + l,i + 2 and so on up the 1 in 
position j,j that is shifted to position j, i. The determi- 
nants of the two kind of matrices are easy to determine 
and they are given by 



and 



det{E{i,j)) = -1 



det{s{i,j)) = {-iy-\ 



The successive apphcation of the matrices E{i + k,j + k) 
for k = up to k = N defines a matrix 



N-l 



Pn{i,j) = n E{i + k,j + k) 



fe=o 



that exchanges a set of A'' columns at once 



di+N-l 



dj+N-l 



PN{i,j) 



"3 ■ ■ ■ 



Applying the matrix to the left of the matrix the 
corresponding exchange of rows is produced. As a con- 
sequence 



PN{i,j) 



AN{i,j) ■■■ AN{i,i) 



V 



(1) 



where Ai^{i,j) are sub-matrices of dimension N x N 
whose first element is located in the row i and column j 
of the matrix where the An are embedded. The result 
obtained, together with that of the pfaffian of a bipar- 
tite matrix will be useful to reduce some of the pfafRans 
encountered in the main body of the paper. Obviously 
detP/v = (—1)^- Unfortunately, this trick can not be 



applied when the number of columns to be "exchanged" 
is not the same. In such case, we have to consider the 
more general shift operation S{i,j). To see how it works 
let us consider a matrix A were there are three groups 
of columns denoted by L, C, and R such that the first 
group goes from column i to columns i + — 1 (i.e. Ni^ 
columns), the group C goes from column i + Nl up to 
columns i + Nl + Nc — 1 (i.e. Nc columns) and finally 
the group R from column i + + Nc up to column 
i + Nl + Nc + Nr — 1 (i.e. Nr columns) . Schematically, 
the columns of the matrix A could be represented as 

A = {...\L\C\R\...) 

Now we want to exchange the group L with the group R; 
to do so the group R of columns is moved to the position 
occupied by L using the product of matrices 



P, 



(1) 



LR 



Nn-1 

H S{i + k,i + k + NL + Nc) 

k=0 



givmg 



{...\l\c\r\...)pI^I^{...\r\l\c\...) 



Now, the group of columns C is moved to the position of 
the group L by means of the following product of "shift" 
matrices 

Nc-l 

Plc= n S{i + k + NR,i + k + NR + NL) 

fe=0 

Using this matrix we obtain 

(. . . \L\C\R\ . . .) pWp^c = (. . . \R\L\C\ . . .) Plc 

= {...\R\C\L\...) 

The matrix exchanging the set of columns L with the 
set R will be denoted Plr = P^^Plc- By applying Pj^ 
to the left of the matrix the set of rows L and R are 
exchanged. As a consequence 



Plniid) 



\ 



■ ANLXNLihi) ■■■ ANLxNniiJ) 

■ ANRxNL{j,i) ■■■ AnrxMrUJ) 



PLR{i,j) = 



•• AnrxMrUJ) ■■■ ANnxNL{j,i) ■■■ 



ANLXNuiiJ) 



i-NLXNL 



V 



where Aj^j^xNLih *) is a sub-matrix with Nl rows and Nl 
columns starting at row i and column i, A^^xNaih j) is 
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a sub-matrix with iV^ rows and Nji columns starting at 
row i and column j — i + Nj^ + Nc, and so on. This 
result generalizes the one of Eq ([T]) for sub-matrices of 
different sizes. The impact of such exchange of rows 
and columns in the pfaffian is the product of the de- 
terminants of the S matrices involved, i.e. det(PLfl) = 
det(P^'j) det (Plc). With det(Pl'^) = {-i)(Ni-+^o)Nn 
and det(Pi,c) = {-1)^lNc ^he total phase is det(PLfl) = 
(^_l-)N^Nn+NcNn+N^Nc that reduces to (-1)^ when 
Nl — iVfl = iV, a result that is independent of A^c"- 



II. THE PFAFFIAN OF A BIPARTITE MATRIX 

To derive some of the results obtained in this paper, 
it is often required to compute the pfafhan of a bipartite 
skew-symmetric matrix with the general structure 



S = 



M Q 

-Q^ N 



(1) 



where M and N are skew-symmetric matrices and Q is 
a general rectangular matrix. Using Aitken's formula it 
is possible to diagonalize the bipartite matrix using a 
congruence transformation 




M 

N + Q^M-^Q 




or the equivalent expression 

I o \ _ 
n-^qt I J ~ 

V n) 

to be used if does not exist. This is a very con- 

venient block diagonalization formula as it involves con- 
gruence transformations that allow to use the property 
pf(P^PP) = det(P)pf(P) of the pfaffian to obtain the 
next two identities 

pi{S) = pf(M)pf(iV + Q^Af-iQ) (2) 
= pf(Af + QiV-iQ^)pf(iV) (3) 



III. THE OVERLAP MATRIX 

The overlap matrix M has the bipartite structure 



T+ 



T 

1(1) 



where the rectangular matrix T, with matrix elements 
Tfcfc' = {afc(0),a+(l)} = (0|afe(0)a+(l)|0), is the matrix 



of the overlaps between the elements of the two basis 
considered. As in the body of the paper, Af(o) and Af(i) 
denote the dimensions of the each of the basis and it is 
assumed for definiteness that TVj-Q-) > iV(i-).The matrices 
1(0) E^iid stand for the identity matrices of dimen- 
sions iV(o) E^nd "^(1)' respectively. For the developments 
considered in this paper, the analysis of the spectral de- 
composition of the overlap matrix is required in order 
to handle properly the occurrence of very small or zero 
eigenvalues of the overlap. The analysis is based on the 
Singular Value Decomposition (SVD) [2^ of the matrix 
T 



T = E+AF 



(1) 



where E and F are square unitary matrices of dimen- 
sions 7V(o) X A''(o) and 7V(i) x A''(i) respectively and A is 
a rectangular matrix of dimension A'^^q) x iV(x) with the 
"diagonal structure" 



A 







where A is a real and positive square diagonal matrix 
with dimension iV(i) x -/V(i). It is convenient to introduce 
a rectangular "identity matrix" I(oi) of dimension A'^(o) ^ 
N(^i^ with a structure similar to the matrix A, namely 



and such that A = I(oi)A. This matrix also has the 
property I(oi)I(oi) = I(i)- 

Using the SVD of T defined in Eq. ([T]) we can finally 
write 



V ^+ 1(1) y 



where 



D = 



E+ 
F+ 



The matrix in the middle can be easily brought to diag- 
onal form 

£)(0)+ ( 1(0) A \ ^(0) ^ 

V A+ 1(1) ; 

/ 21(0) +I(oi)(A- 1(1) \ 
V I(i)-Ay 

by means of the P''-^' transformation 



a/2 ( 1(^01) 1(1) 



£)(0) = J_ I ^(0) "^(01) 
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Finally, introducing the matrix D = DD^^^ we obtain 
the complete diagonalization of the overlap matrix 



21 



(0) 



1(01) 





Id))! 



(01) 







Hi) 



A 



... (2) 

The semi-positive character of the matrix M implies that 
all its eigenvalues are positive or zero and therefore (lower 
block) the elements of the diagonal matrix A cannot ex- 
ceed one. In the case of bases with equal dimensions 



(0) 



the above result becomes 




where the diagonal matrices n± = I± A have been intro- 
duced. As the matrix A is also positive semidefinite it is 
clear that the zero norm eigenvalues of M are associated 
to values of A equal to one. 

It is helpful to consider the special case when T is an 
unitary matrix. It corresponds to the common situation 



where the two bases are connected through a unitary 
transformation as, for instance, the ones associated to 
symmetry operations like rotations in real space, etc act- 
ing on closed basis. For unitary T, the SVD factors can 
be chosen as = T, A = I and F = I in the formulas 
above so that 



D 




(4) 



and 



N=D 



21 




(5) 



with iV-fold degenerate eigenvalues 2 and 0. 

To finish this appendix, just mention that a similar 
(but less general) treatment of the norm overlap was con- 
sidered in Ref [29| but without resorting to the powerful 
concept of the SVD. 
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